Stochastic Resonance in the Fermi-Pasta-Ulam Chain 
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We consider a damped /3-Fermi-Pasta-Ulam chain, driven at one boundary subjected to stochastic 
noise. It is shown that, for a fixed driving amplitude and frequency, increasing the noise intensity, 
the system's energy resonantly responds to the modulating frequency of the forcing signal. Multiple 
peaks appear in the signal to noise ratio, signalling the phenomenon of stochastic resonance. The 
presence of multiple peaks is explained by the existence of many stable and metastable states that 
are found when solving this boundary value problem for a semi-continuum approximation of the 
model. Stochastic resonance is shown to be generated by transitions between these states. 
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Since its theoretical discovery P, [3] and experimental 
verification [^, 0] stochastic resonance has become one 
of the most spread topics in many different branches of 
physics, biophysics and chemistry (see Ref. Q for a re- 
view). It has a wide variety of appUcations, e.g. in semi- 
conductors, in neuronal systems, electronic and magnetic 
systems and, recently, even in quantum physics. All these 
studies are unified by the common idea to represent a 
bistable system as particle motion in a double well poten- 
tial. According to this simplification the system begins 
to jump from one well to the other, following the periodic 
forcing, when a resonant noise level is attained. These 
concepts have wide applications to optical systems. In 
the ring laser, bistability is due to the left-right propaga- 
tion symmetry of the system, while in absorptive optical 
media it is a consequence of the existence of different 
output powers corresponding to the same input pump 
intensity Q. This latter situation has some similarities 
with the system we discuss in this Letter. 

This approach cannot be straightforwardly applied to 
spatially extended systems, where different stable and 
metastable excited regimes can coexist, expecially when 
the forcing is applied locally. As discussed in Q , the pres- 
ence of metastable states does not hinder stochastic reso- 
nance. In this Letter, we show that stochastic resonance 
can be realized for the well-known Fermi-Pasta-Ulam 
(FPU) chain Q of anharmonic oscillators. Up to now 
stochastic resonance in coupled one dimensional systems 
has been considered only when each element of the chain 
displays bistable properties d, In our case, bistabil- 
ity characterizes the whole chain, being originated by the 
coexistence of different stationary regimes corresponding 
to a single driving amplitude. This coexistence is also 
at the basis of the nonlinear bistability effect [ll|, 111, 
which appears when the chain is harmonically driven at 
one end with out-band frequency and is due to presence 
of breather like excitations in the system. 

The equations of motion for the /3-FPU chain (inter- 
particle potential with quadratic and quartic terms) are 
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FIG. 1: Typical signatures of stochastic resonance. In the 
top-left graph we display the power spectrum of the energy 
signal of the FPU chain for the two noise strengths D = 0.03 
and D = 0.09 at which the signal to noise ratio (SNR) shows 
the marked peaks reported in the top-right graph. The forcing 
signal, with carrier frequency Q — 2.05 and modulation W = 
0.003, is plotted in the central graph. The time evolution of 
the energy for the noise intensity corresponding to the first 
peak of the SNR plot is shown in the bottom graph; it displays 
the characteristic bistable features of stochastic resonance. 



Un = Un+l + Un-1 - 2m„ + (u„+l - M„) 

-f (u„_i - Un)^ - (7Un + S,n{t) , (1) 

where u„ stands for the displacement of the n-th unit 
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mass of the chain in dimensionless units (n = 1, 2 ... , N); 
a is the damping parameter (fixed throughout the paper 
to CT = 0.04) and £,n{t) is a zero average Gaussian white 
noise appUed independently to each oscillator with the 
following autocorrelation function 



{U{t)Cnm = 2D5{t)5r, 



(2) 



We force the system at the left end as follows 

un{t) = Acos(m)cos(m), (3) 



while the right boundary condition is free 

Ujv(t) = UN+l{t) . 



(4) 



As commented below, the value of N should not be too 
large, in all simulations presented here = 6. 

In Fig. 1 we display several characteristic features of 
stochastic resonance in the FPU chain. Let us first con- 
centrate the attention on the top-right plot, where we 
show the signal to noise ratio in decibels. Those who 
are familiar with stochastic resonance may be surprised 
to observe multiple peaks in this plot while, usually, a 
bell shaped curve with a single maximum is found. We 
will explain below this observation, providing an analyt- 
ical argument based on the continuum limit of the FPU 
model. The presence of multiple peaks derives from the 
existence of multiple stable and metastable states for a 
given forcing amplitute A (the equivalent of having multi- 
ple wells in the potential). The forcing signal shown 
in the central graph in Fig. 1, has a carrier frequency 
il which is chosen slightly above the upper edge of the 
linear spectrum fl > ujq ^ 2 in order to excite breather 
type localisations, and is modulated with a much smaller 
frequency W <^ fl. This latter frequency plays the role 
of the forcing frequency appearing in standard stochastic 
resonance. There is a wide range of values of both A and 
W for which we can observe stochastic resonance in FPU, 
given that we respect the important condition tr S> W. It 
must be however remarked that the amplitude A should 
be in a range where a stable stationary state of the chain, 
discussed below, coexists with a metastable one. When 
this condition is met, we vary the noise strength D until 
we observe the typical oscillations of the energy of the 
system shown in the lower graph of Fig. 1. The power 
Fourier spectrum of the energy signal shows a sharp peak 
at the driving frequency W when the noise intensity is 
at resonance (top-left graph in Fig. 1). 

We develop below our analytical approach. Let us be- 
gin by solving the boundary value problem for the FPU 
chain in the absence of damping and noise. It is well 
known that the FPU chain is characterized by the fol- 
lowing linear dispersion relation oj{k) = ^^2(1 — cos fc), 
hence the spectrum has band edge frequency ujq = 2. In 
order to derive stationary weakly nonlinear solutions (see 
e.g. Ref. [l^) of Eq. ([T|), we seek for solutions of the 



following standard form [1^, 17 1: 

(-1) 



Unit) 



(5) 



where c.c. denotes the complex conjugated term and 
(j){n, t) varies slowly with respect to both its arguments 
n and t. These solutions should be considered as mod- 
ulations of upper band edge oscillations. The boundary 
value problem ^ can be rewritten for 4>{n, t) in the form 

0(O,i) = Ae*(^-2)* 0(iV + l/2,t)==O, (6) 

where the second condition derives from (j){N -f- 1) — 
—(t){N). Moreover, the condition Sui = fl — 2 1 must 
be satisfied for the function 0(n, t) to be slowly varying 
with time. 

Substituting then ^ into Eqs. ([T]), and assuming 
(^(n, t) to be a continuous function of its variables, one 
gets, in the weakly nonlinear limit and neglecting higher 
order derivatives in n and t, the following nonlinear 
Schrodinger equation 3, 17 1: 



dt 



4 dn"^ 



0. 



(7) 



We now assume that the system synchronizes with the 
boundary ^ and define a new parameter B 



(t>{n,t) 



AtSuj 



(p{r 



B = 



dip{n) 



dn 



(8) 



n=N +1/2 



(where Lp{n) is a real function of n) and, as a result, we 
get the following relation 

|^)' = s2+4(f^-2)^2_g^4 

and its solution in terms of Jacobi elliptic functions 
1^ , isj . which, after substitution into Eqs. and (O 



gives the following approximate solution for the relative 
displacements of the oscillators 



Un = (-l)"Qcos(r2t) cn 



2^[N+^-n] -IC{k), k 



(10) 

where /C(fc) is the complete elliptic integral of the first 
kind with a modulus k and all the constant are defined 
via the single free parameter B as follows 
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-B' 



Slu + ^ 



3 47(7 — Sllj) 

. . 

The free parameter B is defined by the condition of 
adaptation to the boundary ([3]), leading to the following 
consistency relation 



A^Qcn [2y/^ {N + 1/2) - IC{k), k] 



(12) 
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FIG. 2: In the top-left plot the energy of the approximate 
solutions (|10|l of the FPU chain obtained using our semi- 
continuum approach is displayed as a function of the driving 
amplitude A. The intersections of the vertical solid line at 
A — 0.09 with the different branches of the curve correspond 
to multiple solutions of the consistency relation (|12|l . Type 
1 {El ~ 0.1) and 2 {E2 — 0.6) solutions, which turn out to 
be stable, are compared with numerical simulations (points) 
in the top-right and bottom-left panels. A type 3 {E3 ~ 1.6) 
solution, which is metastable (see text), is represented in the 
bottom-right panel. Ai = 0.18, A2 = 0.31 and ^3 = 0.44 
are threshold values for the excitation of the different solu- 
tions: Ai is the driving amplitude for which the transition 
from a type 1 stationary state to a type 2 stationary state 
occurs, while A2 marks the transition from type 2 to type 3 
(metastable) solutions, and so on. In the figure, we also dis- 
play a feature which is present because of the damping and 
is not taken into account by our approximate solutions: an 
additional lower threshold L is present (vertical dashed line), 
below which all the solutions collapse to the type 1 lowest 
energy state. 



From this equation one gets all compatible values of B 
for a fixed driving amplitude A. It turns out that there 
are multiple values of B that solve Eq. (H^), which in 
turn implies the existence of multiple compatible pat- 
terns of the type ((TO)) . which have correspondingly dif- 
ferent energies. In Fig. 2 wc display these energies as 
a function of the driving amplitude A. Three branches 
of solutions are present. For instance, considering the 
forcing amplitude A = 0.09, five possible energy values 
are found: among these two are unstable, two are sta- 
ble and one is metastable, as described in the caption 
of Fig. 2. For the type 1 and 2 solutions in Fig. 2, 
one has an excellent correspondence between the ana- 
lytical expression (jlOp and the direct numerical simula- 
tions of the FPU chain. For type 3 solutions the approx- 
imation is not as good and, besides that, the solution is 
metastable: in numerical simulations it destabilizes and 
restabilizes during time evolution. We do not display in 
Fig. 2 solutions corresponding to higher energies, be- 
cause the semi-continuouum approximation fails and one 




FIG. 3: Time evolution of the energy of the FPU chain for a 
constant driving amplitude A — 0.18 and a linearly varying 
noise intensity. Transition from a type 1 to a type 2 state is 
shown to occur when the noise intensity reaches the threshold 
value D — 0.021, while the transition to a type 3 states takes 
place at _D = 0.055. In this simulation, damping is present 
and the driving is not modulated, W = 0. 

should use truncated wave approximations [l^ . 

Summarizing, the system is characterized by several 
threshold amplitudes. When the driving amplitude ex- 
ceeds the first threshold Ai, the type 1 stable state jumps 
into the type 2 stable state. By further increasing the 
driving up to the second threshold A2, the transition to 
the type 3 metastable state occurs. In the presence of 
damping, when reducing the amplitude of the driving, 
the system does not remain in the higher energy states. 
Indeed, a lower threshold amplitude appears, L = 0.03, 
below which the system goes back to the lowest energy 
stable state. Afterwards, when we will consider stochas- 
tic resonance, we will modulate the driving signal in such 
a way that its amplitude remains above L and below Ai 
(typically A G [0.10,0.14]), in order to realize stochastic 
transitions between all states. Furthermore, we note that 
for a fixed driving amplitude A the stable and metastable 
states are characterized by sufficiently far separated en- 
ergies. The presence of these "energy levels" can be mon- 
itored by keeping the driving amplitude constant and 
changing the noise intensity, as shown in Fig. 3. It is 
clearly seen that for the fixed boundary driving ampli- 
tude A = 0.18 the transition between type 1 and type 
2 stable states occurs at the noise intensity D — 0.021, 
while at D = 0.055 the system goes from a type 2 sta- 
ble state to a type 3 metastable state. One should also 
observe that the averaged energies of each state well cor- 
respond to the ones derived from the approximate ana- 
lytical solutions ((To)) . 

Now we are ready to discuss the main result of this 
paper. Simulating the FPU model ([T|) with the driv- 
ing signal seed frequency f2 = 2.05 and modulating fre- 
quency W = 0.003, such that the driving signal ampli- 
tude A varies in the range [0.10,0.14], in the presence 
of both damping and noise, we monitor the time evolu- 
tion of the system's energy for increasing noise intensities 
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FIG. 4: Time evolution of the signal (top) and system's en- 
ergy (central and bottom) at the noise intensities D — 0.05 
(corresponding to the point between first two peaks of the 
SNR in Fig. 1) and D = 0.09 (second peak of the SNR), 
respectively. 



D, and then analyse the properties of its power spec- 
trum P{iLj). The signal to noise ratio (SNR) is defined as 
SNR = 10logi„{P{ujpeak)/N{ujpeak)), where ujpeak is the 
frequency at which the power spectrum displays a clear 
peak (see the top- left plots in Fig. 1) and N{ujpeak) is 
the power of the background noise at the same frequency. 
For low noise intensities there is no transition between the 
different states and the system's energy oscillates around 
El. For the noise intensity value D = 0.03 (bottom graph 
in Fig. 1) one observes an oscillation between the state 
with energy Ei ~ 0.1 and that with energy E2 — 0.6. In- 
creasing further noise intensity, the SNR drops and the 
energy fluctuates around the value E2 = 0.6 (middle plot 
in Fig. 4 for D = 0.05)). Finally, for the noise intensity 
D = 0.09, which corresponds to the second peak of the 
SNR in Fig. 1, we observe noisy transitions between the 
metastable state of energy E3 = 1.6 and two stable states 
with energies Ei and E2 (see the bottom graph of Fig. 4). 
Wc do not discuss here the smooth third peak observed 
in Fig. 1, which should be connected to transitions with 
higher energy states, of which we do not possess a clear 
analytical characterization. 

As remarked above, we have performed all our numer- 
ical analysis for short FPU chains (throughout the paper 
N = 6), but the observation of this phenomenology is 
not restricted to this case. Indeed, it should be men- 
tioned that, if in the equations for the /3-FPU chain ^ 
we would have rescaled the strength of the nonlinear cou- 
pling to /, the semi-continuum approach developed here 
would have produced exactly the same solutions PTI)) . 
provided that the number of particles N would have been 
increased maintaining the ratio \/f /N constant. 

Concluding, we have considered a damped /3-FPU 
chain forced at one boundary with a modulated signal 



under the action of an increasing noise level. For specific 
values of the noise intensity, the power spectrum of the 
system's energy displays sharp peaks at the modulating 
frequency. Correspondingly, also the signal to noise ra- 
tio shows pronounced peaks as a function of noise inten- 
sity. This is a clear observation of stochastic resonance 
in an extended system. In contrast to all previous stud- 
ies, multiple peaks of the SNR are observed, instead of a 
single broad bump. In order to explain this feature, we 
have developed a semi-continuum analytical approach, 
which allows to point out the presence of multiple stable 
and metastable solutions for a given forcing amplitude. 
Stochastic resonance is shown to be generated by transi- 
tion between these states. These studies could be readily 
extended to realistic physical systems with local forcing. 
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